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3 He and the triton are studied as three-body bound states in the effective field 
theory without pions. We study 3 He using the set of integral equations developed 
by Kok et al. which includes the full off-shell T-matrix for the Coulomb interaction 
between the protons. To leading order, the theory contains: two-body contact in- 
teractions whose renormalized strengths are set by the NN scattering lengths, the 
Coulomb potential, and a three-body contact interaction. We solve the three cou- 
pled integral equations with a sharp momentum cutoff, A, and find that a three-body 
interaction is required in 3 He at leading order, as in the triton. It also exhibits the 
same limit-cycle behavior as a function of A, showing that the Efimov effect remains 
in the presence of the Coulomb interaction. We also obtain the difference between 
the strengths of the three-body forces in 3 He and the triton. 



INTRODUCTION 



Since Weinberg first proposed applying the ideas of effective field theory (EFT) to nuclear 
forces [H, much effort has gone into this approach. (For reviews, see Refs. [2-4|.) Although 
there is still some debate about how best to implement it at energies where pion-exchange 
forces are resolved, the picture is clearer at lower energies. Here few-nucleon systems can 
be described by a "pionless" EFT based on two- and three-body contact interactions js-8]. 
The resulting expansion of the two-body force is just that of the effective-range expansion 
(9)], but the EFT framework makes it possible to extend this consistently to other effective 
operators and three-body forces. 

This theory has been applied extensively to two-body systems, where it has been extended 
to include the effects of the Coulomb interaction on proton-proton scattering 10Nl4|. In 



that system, it corresponds to a distorted-wave or "modified" version of the effective-range 
expansion [9]. 



Three-body systems have also been studied in the pionless EFT 0, Q] . In these, exchange 
of one particle between an interacting pair and the third particle leads to a long-range force 
that can be either attractive or repulsive depending on the overall symmetry of the wave 
function. Attractive channels, such as the triton, display the Efimov effect 15] in the limit 
of infinite two-body scattering length. This consists of an infinite tower of bound states with 
energies in a constant ratio. It results from a discrete scale invariance in these systems and 
corresponds to a limit cycle of the renormalisation group (RG) for the contact three-body 
force [16 1 . 

The 3 He nucleus is of particular interest since it can potentially provide access to prop- 
erties of the neutron that require a polarised target. It should be possible to describe its 
low-energy properties within the pionless EFT, extended to include the Coulomb interac- 
tion. However the only previous application of the theory to the pd system is in the work 
of Rupak and Kong, who treated Coulomb effects perturbatively and only considered the 



spin-3/2 channel [17 ]. 



A nonperturbative treatment may not be essential for the 3 He ground state, given the 
typical momenta in the wave function. However it will be necessary if the EFT is to be 
extended to describe pd scattering near threshold. Also, it is needed to answer the question 
of whether the Efimov effect survives in the presence of the Coulomb interaction. Arguments 
based on the degrees of the singularities in the potentials suggest it should and this has 



been checked by Hammer and Higa for a simpler two-body model 18j. However it has not 
previously been confirmed in a three-body system. 

Like most other EFT studies of the triton and other three-body systems ]3], our treatment 
is formulated in terms of a set of integral equations. However a quite different approach, 
based on the resonating-group method, was recently applied to three- and four-nucleon 



systems by Kirscher et al. 19]. While this also treats the Coulomb interaction to all orders, 
limitations on the range of cutoffs mean that it is not able to test whether the Efimov effect 
occurs. 

Here we report the results of a study of 3 He in the pionless EFT with the Coulomb 
interaction treated nonperturbatively. Our approach is based on that used by Kok et al. 
{20! to study 3 He with separable potentials and the full off-shell Coulomb T-matrix. The 



resulting set o 



integral equations is very similar to that developed by Skornyakov and ter- 



Martirosian 21] which have more recently been used in EFT studies of the triton and other 
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three-body systems [3]. 

Our results show that the Efimov effect does indeed occur in the presence of the Coulomb 
interaction. We also determine the strength of the three-body force needed to reproduce 
observed 3 He binding energy. The difference between this and the corresponding force needed 
for the triton is of the expected size for an electromagnetic effect. 



II. LAGRANGIAN 



In the pionless EFT, the strong forces between nucleons are described by two- and three- 
body contact interactions. In the present context it is convenient to represent the two-body 
interactions in terms of dibaryon fields. The resulting effective Lagrangian for the three- 
nucleon system can be written as jsl, 22^ 



C = C 



A' 



C s + Ct + £3 , 



(1) 



where Cn is the standard one-nucleon Lagrangian in the heavy-baryon formalism, 

£ v = {iv ■ D + \(v ■ D) 2 -D 2 ]\N. (2) 



Here is a velocity vector satisfying a condition v 2 = 1, is the covariant derivative, and 
m^v is the nucleon mass. The terms C s and C t are dibaryon effective Lagrangian for spin 
singlet and triplet parts, respectively, and these read 



C s = s{{iv- D + A a ( a ) } s a - y s \ s{ 
C t = t\ {iv D + At} ti-y t (t\ 



N T p£ So) N 

+ h.c. J 



h.c. 



N T P? Sl) N 



(3) 
(4) 



where s a and £, are the corresponding dibaryon fields. The strengths of the corresponding 
two-body interactions depend on A s ( ffl ) and A f , the mass differences between the dibaryons 
and two nucleons, and y s and yt, the coupling constants for the dibaryon- nucleon- nucleon 
vertices. The projection operators for the two-nucleon 1 5'o and 3 Si states are 



P( 3S1 ) - 1 T a n 



(5) 



where r a and (jj are Pauli matrices for isospin and spin, respectively. 

The present work is based on the leading-order (LO) terms in the expansion of this 
EFT of Q. Here Q stands for any of the low-energy scales in the problem: the relative 
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momenta, 1/a where a is any of the NN scattering lengths, and the inverse of the Bohr 
radius k = a EM m N /2. More precisely, we keep all terms of order Q^ 1 , noting that the 
strong attraction in the NN S waves enhances the wave functions near the ori gin and hence 



promotes contact interactions proportional to 1/a or k to this order p, [D, |23|, \24 \. 

The Coulomb interaction leads to a well-known logarithmic divergence in the NN loop 
integrals 1CH12|. This is renormalised by a LO counterterm that contributes to the pp 
scattering length. Given that we need to include isospin breaking in this channel, we have 
also taken account of charge-independence breaking between the 1 Sq np and nn scattering 
lengths, even though this is formally of higher order. This allows us to present the general 
structures of the three-body equations, which will be needed for extensions of this work to 
higher orders. We have not included the neutron-proton mass difference which is not only 
of higher order but is smaller than might be expected as a result of cancellations between 
electromagnetic and strong contributions 25|. To the order we work, observables depend 
only on the combinations A s ^/y 2 and so we have chosen to leave the couplings y s isospin- 
invariant. 

The three-body force is expressed as a dibaryon-nucleon contact interaction. It is given 
by the effective Lagrangian, 

+ y 2 s N\r-s)\r-s)N) ) (6) 

where H(A) is the coupling constant, which runs with the scale A of the cutoff we impose 
on the coupled integral equations. Note that this interaction contributes only to three-body 
channels with total spin 1/2 and total isospin 1/2. 

The building blocks needed to construct the three-body integral equations from the La- 
grangian are as follows. In the one-body sector, we have the non-relativistic nucleon propa- 
gator, 

iS N (p) = a -. (7) 



In the two-body sector, we need the dressed dibaryon propagators with and without the 
Coulomb interaction. The propagators for the np and nn channels, which have no Coulomb 



FIG. 1: Diagrams for the dressed dibaryon propagator, denoted by the double line with filled circle. 
Single lines represent nucleon propagators; double lines undressed dibaryon propagators. 



interaction, are represented by the diagrams in Fig. [T] and are given by 

iD Stt E - -2— = - 2 . 8 

V 4m N J m N y Sit _^_^_ + ^3 q2 _ mNE _ %e + te 

The two-nucleon loop diagrams here have been dimensionally regularised using the PDS 
scheme with the subtraction scale //. The dependence on this can be removed by renormal- 
ising the constants 47rA Si t/ (mjvyft) using 

47rA s („p) 1 47rA s(nn) 1 47rA t 

/i + V = , H+ ~t 1 = , yu + 2=7, (9) 

m N yj a np m N yj a nn rn N yf 

where a np and a nn are the scattering lengths for the spin-singlet np and nn channels, respec- 



tively, and 7 is related to the deuteron binding energy B 2 through 7 = y/m N B 2 . Working 
to leading order, we neglect corrections from the effective ranges. 






FIG. 2: Diagrams for dressed dibaryon propagator for the pp channel, denoted by the double line 
with shaded circle. The shaded ovals here represent the two-nucleon Green's function dressed with 
the Coulomb interaction. 

In the pp channel, the Coulomb interaction dresses the two-nucleon Green's function in 
the bubble diagram for the dressed dibaryon propagator, as shown in Fig. |2j The resulting 
propagator is given by 

3q 2 \ 4vr -i 



where 



iD s(pp) E - — — = —— — ^ , (10) 

4m N J m N y s _ _ 2kH (k/ P ') 



H{rf) = V(^) + ^~ - Hirj) , (11) 
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ip is the logarithmic derivative of the T function, and — ip' = y |g 2 — ra^E — ie. The 
renormalised constant A^ pp ^ for the pp channel cancels both the linear and logarithmic 
divergences of the loop diagram. It is related to the Coulomb scattering length ac by 



1 

a c 



4ttA 



1-C f 



In 



4k 



(12) 



where Ce = 0.577215 • • • is Euler's constant. Note that the logarithmic divergence means 
that it is impossible to make a model-independent decomposition of 1/ac into strong and 



electromagnetic contributions 10Ml2| . This implies that within this EFT there is no unam- 



biguous way to separate Coulomb from other isospin-breaking effects. 

The final building block for the integral equations is the off-shell Coulomb T-matrix. A 



convenient form for it is the integral representation 



26j for negative energies, k 2 /m N < 0: 



(p'\T c (k 2 /m N )\p) 
where rj = n/k and 



IP " P 



4 "'/ -«)•(,. -i; 



(13) 



X 



1 + 



(p' 2 - k 2 )( P 2 - k 2 
jfc2(p'-p)2 



(14) 



III. INTEGRAL EQUATIONS 

Now we construct the integral equations for scattering of a third nucleon off a deuteron, 
concentrating on the channels where the third nucleon is in an S wave, and the total spin and 
isospin of the three particles are both 1/2. The negative-energy solutions of these describe 
the bound triton and 3 He states. 

We work the center of mass frame for the Nd scattering and use the notation of Bedaque 
et al. jsj]. In the isospin-symmetric case, this process can be described by an amplitude 
a(p, k) for scattering into states with a 3 Si dibaryon, and b(p, k) for scattering into those 
with a 1 Sq dibaryon. Here k denotes the initial (on-shell) relative momentum of the deuteron 
and the third nucleon, and p the final (off-shell) momentum. 

The amplitudes satisfy coupled integral equations which can be represented diagrammati- 
cally as in Fig. [31 The kernels of these equations consist of two terms: one-nucleon-exchange, 
which provides a long-range force between the dibaryon and the third nucleon, and the three- 
body contact interaction. 
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In the limit of infinite two-body scattering lengths (and zero deuteron binding energy) 
the long-range force is scale- independent. It is this feature that leads to the Efimov effect 
in channels where the force is attractive. 




FIG. 3: Diagrams corresponding to the coupled integral equations for nucleon-dibaryon scattering. 
Thin lines denote nucleons, double and thick lines 3 Si and 1 So dibaryons. Those with filled circles 
are dressed dibaryon propagators, as in Fig. [TJ Ovals with stripes denote the off-shell amplitude 
a{p, k) with initial and final spin-triplet dibaryons, and those with cross stripes the amplitude 
b(p, k) with initial spin-triplet and final spin-singlet dibaryon fields. 

When we include the Coulomb interaction, isospin symmetry is broken and we must 
introduce separate amplitudes b+(p,k) and bo(p,k) for pd scattering into states with 1 S'o 
pp and np dibaryons, respectively. These and a(p, k) satisfy a set of three coupled integral 
equations. For completeness, we also present the corresponding set of three equations that 
arise when we allow for isospin-breaking effects in the nd case, involving the amplitude 
6_ (p, k) for states with nn dibaryons. As already mentioned, these are not essential for the 
present LO calculations but will be needed for extensions to higher orders. 

A. Integral equations for the nd channel 

Before presenting the equations for pd scattering, we first look at the simpler equations 
for nd scattering. Allowing for isospin breaking in the scattering lengths, the three integral 
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equations for the S = 1/2, T = 1/2 channel are 
a(p,k) = K^(p,k;E)+2^ 



A 2 



where 



+ 



dq q 2 



dq q 2 



dq q 2 





a{q,k) 


A 2 J 


-7 + 


fa - 


mE 


2 if (A)" 




bo(q,k) 




3 A 2 


i_ 




— mE 


2 ff (A)" 




b-(q,k) 




3 A 2 


1_ 


+ \/h 2 


— mE 



b (p,k) = 3K^ a \p,k;E)+2 1 ^ 

-A 



+ 



7T 



+ 



7T 



+ 



7T 



c/g q 



dq q 2 



dq q 2 



3K^(p,q;E)+2^- 
-K^(p,q;E) + 2 -^ 
K^(p,q;E) + 2 -^ 



a(q,k) 



-7 + ^\q 2 - mE 
b (q,k) 



|g 2 — mE 



b_(q,k) 



6_(p,fc) = 3tfW(p,fc;£)+2^ 



+ 



7T 



+ 



7T 



2 

+ - 

7T 



dq q 2 



dq q 2 



dq q 2 



3K^(p,q;E)+2^- 
'^\p,q;E) + 2 -^ 



a(q,k) 



-7 + y |g 2 - 
6o(g,A;) 



- + \/i<? 2 - m ^ 



2 ff (A) 

3 "A 2 " 



b_(q,k) 



— + \ \q 2 -mE 



^(a)/ _g,N = 1 \ n ( p2 + + 2 PQ- m NE \ 

2pq \p 2 + q 2 — 2pq — m^E J 



(15) 



(16) 



(17) 



(18) 



If we set a nn = and b = b_ = b, these reduce to the isospin-symmetric forms of the 
equations in Ref. [8]. 

The integrals over the relative momentum q are all cut off at q = A. The resulting 
dependence on A can be cancelled by the three-body force, whose strength is if (A)/A 2 . 
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B. Integral equations for the pd channel 

We now turn to the corresponding equations for pd scattering. These differ from the ones 
in the previous subsection by having b + instead of b- and, more importantly, adding the 
Coulomb interaction between the two protons. They can be obtained from the equations 
developed by the Groningen group [20| , by omitting the separable form factors and instead 
imposing a sharp cutoff on the relative momenta. 

The Coulomb interaction leads to additional long-range terms in the kernels of these 
equations. The diagrams for the full long-range kernel are shown in Fig. HI The one-nucleon 
exchange term of diagram (a) is supplemented by diagram (b) if a proton is exchanged, 
and by (c) for neutron exchange. Finally, diagram (d) represents the Coulomb interaction 
between an np dibaryon and a proton. 




(a) (b) (c) (d) 



FIG. 4: Diagrams corresponding to the long-range terms in the kernel for pd scattering. Here, the 
shaded oval denotes the off-shell Coulomb T-matrix (since contributions where the two protons do 
not interact are already contained in (a)). Diagram (b) has been shown for the case of an initial 
np dibaryon and a final pp one. The corresponding diagram for the inverse process is the mirror 
image of this. 
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The coupled integral equations can be written 

a(p, k) = K^(p, k; E) + K®(p, k; E) + 2K^(p, k; E) + 2 



A 2 



dqq 2 



(p, q; E) + (p, k- E) + 2K® (p, k; E) + 2 

A^ 



a(q,k) 



-7 + i/ fg 2 - mE 



+ - [ A dq q 2 f Jf W (p, g; J5) + iT (c) (p, fc; E) + \ ^ 



b (q,k) 



+ - / dqq 1 

71 



K^\p,k;E) + KS(p,k;E) + 2 -^ 



— - + \ \q 2 -mE 
b+(q,k) 



i - 2k#(k/p') 



(19) 



fcoM = 3i^(p,£;;£) + 3K( c )(p,A;;£) + 2^ 

-A 



+ - / dq q 2 (p, q;E) + 3K^ (p, q; E) + 2 ^ 

^ Jo A 



a(q, k) 



<ig g 2 



ff« (p, g; E) + ffM (p, g; E) + 2K® (p, g; £) - | ^ 



b (q,k) 



\q 2 — mE 



+ 



71 



dqq 2 



K^(p,q;E)+K^(p,q;E)+ 2 -^ 



b+(q,k) 



__L _ 2kH{k/p>) 



(20) 



6+(p,A;) = 3#W(p,A ; ;£)+3tf£ ) (p,A;;^ + 2^ 



+ 



dgg 2 



dqq 2 



3#W(p, q; E) + 3J$(p, g; £) + 2 ^ 



a(g, fc) 



2 #(A) 



3 A 2 



- -7 + y fg 2 - wi-E 

&o(g,fc) 

+ \/h 2 - mE 



71 



+ - / dqq 2 



2 #(A) 

3 A 2 



b+{q, k) 



- 2«H(«/j/) 



(21) 



where is the one-nucleon exchange kernel defined in Eq. (jl8jl above, and If^, ^{3 , 
i^ c \ and are the pieces corresponding to the Coulomb diagrams (b), (c), (d) in Fig. HI 
(In and , the subscripts 1 and 3 refer to initial and final dibaryons, 1 denoting np 
and 3 pp.) 
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In diagram (b), the dibaryon- nucleoli- nucleoli vertex projects out only the 5- wave part 
of the Coulomb T-matrix. This considerably simplifies the corresponding terms in the 
equations by reducing the number variables that need to be integrated from six to two. The 
corresponding terms in the kernels are 

<f>fi(p,q,y;E) 



1 



dy 



2 J _i p 2 + q 2 — mE + pqy ' 



1 



dy 



2 J _i p 2 + q 2 — mE + pqy ' 



(22) 
(23) 



where y = p ■ q = cos 9 and 



dtt K ^' 



Ak/3" 



dw/? n 



Here we have defined the momentum variables 



(24) 
(25) 



V 



1 



P+ 2 q 



p 



1 



^q 2 -mE, /3" = J ^p 2 - mE . (26) 



The kernels for diagrams (c) and (d) involve six dimensional integrals but we can analyt- 
ically evaluate two of these (over the azimuthal angles). In the case of K^°\ it is convenient 
to introduce the vectors, 



a = p + q, b = p-q, 



Choosing axes where a lies along the z-direction, these can be written 



(27) 



a = a(0,0, 1), b = &(sin#",0,cos#"), 1 = I (sin0' cos 0', sin 6' sin 0', cos 9') , (28) 



in terms of their magnitudes 



a = \/ p 2 + q 2 + 2pq cos 9 , b = \/p 2 + q 2 — 2pq cos 9 , 



(29) 



and the angles defined by 

2pq sin 9 



sin 9" 



a/ (p 2 - q 2 ) 2 + Ap 2 q 2 sin 2 9 ' 



COS0" 



P 2 ~q 2 



\J (p 2 - q 2 ) 2 + Ap 2 q 2 sin 2 9 



(30) 
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The kernel then becomes 

K (c \p,q-E) = -- I" dl f 1 d(cos9) f 1 d(cos9') f dt {K/P)t K,fi 1 

Wo J-i J-i Jo bl^Jd 2 + 4t(l-t)~ 2 l 2 f3 2 



1 1 



cos 9" cos 6> + l ^d 2 + At(l-t)~ 2 l 2 (3 2 \[Y^y{ 

1 1 



cos 9" cos 9' - ^d 2 + U{l-t)- 2 l 2 P 2 y/T^yl 



(31) 



where 



2/1,2 



sin 9" sin#' 



cos 9" cos 9'±l ^/d 2 + U{l -t)~ 2 l 2 F 



3 

d = I 2 + p 2 + q 2 + pq cos 9 — -al cos 9' — m N E , 
(3 2 = ^(l 2 + a 2 + 2al cos 9') -m N E . (32) 

The last term in the kernel, arising from diagram (d), requires the most care in its numer- 
ical evaluation. This is because the long-range photon exchange gives rise to a logarithmic 
IR singularity. In the present work, where we are interested only in bound states, the finite 
extent of the wave functions helps to regulate this singularity. Calculations of pd scattering 



amplitudes will, however, require different numerical methods. Following Ref. 20(, we split 
up into three pieces, 

(p, q; E) = (p, q; E) + (p, q; E) + K^ T ~ V ^ (p, q- E) , (33) 

where K^ d ' yn> and K^ d ' v "^ denote the contributions the from singular and nonsingular parts 
of one-photon exchange, and K^ d ' T ~ v ^ the remaining (also nonsingular) terms from iterating 
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the Coulomb potential. These have pieces have the forms 



K^ v '\p,q;E) = 



In 



p + q + 2pq 



4pqf3' \p 2 + q 2 — 2pq J ' 
K^ v "\p,q;E) = -« /' d(cos9) 1 

J-i p 2 + <r - 



(34) 



2pqcos6 [ v /p 2 + q 2 -2pqcos9 



x 



arctan 



K^ T - v \p,q;E) = -k 



+ arctan 



J d(cos9) 



—p 2 + 2q 2 — pg cot 6 
2/3'vV + g 2 -2pqcos6 / 

2p 2 — q 2 — pq cos 



2/9 \/p 2 + <? 2 — 2pg cos 
1 

p 2 + g 2 — 2pg cos 



2/3' f ' 



(35) 



x 



dlt 



j d(cos9') 



„ y 1 dt (k/a)^ /a_1 



arctan 



Dv/T^ 2 

— p 2 + 2g 2 pg cos 



a/p 2 + g 2 — 2pg cos 

2p 2 — q 2 — pq cos 



2f3'^p 2 + q 2 -2pqcos6 



+ arctan 



2/3 ^p 2 + q 2 — 2pqcos8 



(36) 



where 



-p 2 - m^E , 



3 /3 

-g 2 - m N E , (3 l = J-l 2 - m N E , 



D = (/ 2 + g 2 + q l cos -m N E)(l 2 + p 2 + pi cos -m N E) 

+4t(l - ty 2 (3f(p 2 + g 2 - 2pgcos0) , 
c = (I 2 + q 2 + ql cos 6-m N E)pl sin sin 0'/-D • 



(37) 



IV. NUMERICAL RESULTS FOR BOUND STATES 

We solve for the three-body bound states by taking the homogeneous parts of the coupled 
integral equations and treating them as nonlinear eigenvalue problems, 

K(E)u = u, (38) 

where u T = (a, b , b±) and the 3x3 matrices of integral operators K(E) can be extracted 
from the expressions in the previous section. 

To convert K(E) into a matrix we use Gauss-Legendre quadrature to replace the momenta 
p and q by discrete variables. We have examined three versions of this discretisation, in terms 
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of variables x defined by: (1) p = mtan(x) where m = 140 MeV, (2) p = x (no change of 
variable), and (3) p = m exp(x) where m = 1 MeV. The variables p and q run from zero 
to the cutoff A, except in the third of these where we also need to impose a cutoff at the 
very small momentum mo- 

We the construct the determinant of the matrix K(E) — 1 and search on E to find the 
energy eigenvalues at which det[K(E) — 1] = 0. For one of these values of E, we then solve 
the eigenvalue problem Eq. (1381) and determine the eigenvector corresponding to the unit 
eigenvalue of the matrix K(E). 

For given values of the two-body input parameters we first set the three-body force 
H(A) to zero and adjust the value of the cutoff until the shallowest bound-state eigenvalue 
reproduces the observed three-body binding energy. Starting from this cutoff, Ao, we then 
determine the values of H(A) needed to reproduce this binding energy for other cutoffs. 



A. Triton 



In this case, we work with the determinant constructed from the homogeneous part of 
the coupled integral equations in Eqs. (fT5l) . (JIB]) , and (fT7|) . The two-body input is provided 
by 



7 = 45.703 MeV, a np = -23.749 ± 0.008 fm [23] , a nn = -18.5 ± 0.4 fm 28] .(39) 
We adjust the three-body force to fit the triton binding energy, 

5 3 ( 3 H) = 8.48182 MeV, (40) 



which corresponds to a momentum scale a = y/m^ B 3 = 89.239 MeV. 

For 16 Gaussian integration points (corresponding to a 48x48 matrix) we find that our 
results have converged and, for the three choices of variable discussed above, they agree to 
6 significant figures. 

For H(A) = 0, we reproduce the experimental binding energy with the cutoff A = 
380.689 MeV. The values of the three-body force needed for other cutoffs are shown in 
Fig. [51 

As we take A — > oo, the cutoff becomes much larger than the physical scales in the 
system (1/a) and H(A) displays the limit-cycle behaviour found in Ref. js]. This reflects 
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FIG. 5: Values of the three-body force H as a function of the cutoff A in MeV. For comparison, 
the dashed curve shows the asymptotic, scale- free form of H(A). 

the presence of the Efimov effect in this three-body system, a tower of deeply bound states 



with energies in a constant ratio 



15| . In three-nucleon systems only the shallowest of these 



states lies within the domain of the EFT and so corresponds to a physical state. 
The asymptotic form of H(A) in this limit is js] 

_ sin[g ln(A/A^ - arctan(l/s )] 
sin[s ln(A/A*) + arctan(l/s )] 

where So = 1.00624 • • • and the scale parameter A* is defined by 

A = A*exp[(l/s )arctan(l/s )] ■ 



(41) 



(42) 



This form is shown by the dashed line in Fig. [5j The numerical results deviate from it only 
for small values of A where the finite scattering lengths can influence the renormalisation of 
the three-body force. 
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B. 3 He channel 

For the 3 He channel, we need to the pp Coulomb scattering length as a two-body input 
parameter, 



a c = -7.8063 ± 0.0026 fm [29j. (43) 
The experimental binding energy of 3 He is 

£ 3 ( 3 He) = 7.71804 MeV , (44) 

corresponding to a scale a = 85.127 MeV 

As already noted, the K^ d,v '^ term contains a logarithmic singularity that needs to be 



dealt with carefully. We handle it with the method outlined in Appendix B in Ref. [20(. In 
the case of a bound state, the other factors in the integral are regular and so the logarithmic 
singularity has a finite integral. We can take the interval around the singular point and, 
treating the other factors as constant, explicitly integrate the logarithm over it. The result 
can be expressed in a similar form to a Gaussian quadrature, 

rb n 

/ dq]n(p j -q) 2 y{q) = ^Wiki(p j -q i ) 2 y(qi)+w' j y{p j ), (45) 

where y(q) is a smooth function and the u>; are the usual weights. The modified weight for 
the singular point is given by 

r<u> 

w j — dq InQo.,- - q) 2 

= -2wj + 2(q j> - pj) ln^ - pj) + 2(p j - q jK ) ln(^ - q jK ) , (46) 

where q 5< = a + YjjZi w k and g 3> = a + J2l=i w k- 

To examine the size of Coulomb effects in the three-nucleon system, we first consider an 
isospin-symmetric three-body force and take H(A ) for the same value A = 380.689 MeV 
as in the triton. We use A4 = 32 Gaussian points for the four integrals needed to evaluate 
the kernels K^ c ' and K^ d \ The integral equations are discretised using up to iV = 32 
points (leading to 96 x 96 matrices). We fit the results for iV = 16 and 32 to the form 
a(N) ~ A + B/N in order to extrapolate to iV = 00. 

Using the first choice for the discretisation variable (p = mtan(x)), we get a = 
84.8084 MeV, or a binding energy of B 3 = 7.66038 MeV. This differs from the triton energy 
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for the same three-body force by 0.82 MeV. It is within 1% of the observed 3 He energy, 
indicating that the isospin-violating three-body force is indeed a higher-order contribution. 
This is as expected, given the absence of any modification of the renormalisation, since the 
isospin-violating term is suppressed by one power of the inverse Bohr radius k relative to 
the LO force. 
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FIG. 6: Evolution with A of the three-body forces H(A) for 3 He (solid line) and the triton (dashed). 



Our difference between the triton and 3 He binding energies is significantly larger than 
the 0.66 ± 0.03 MeV found by Kirscher et al. 19]. However, although working to NLO, 



ength. Their results 
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3l| 



those authors took a charge- independent value for the 1 So scattering 
are consistent with other estimates of the pure Coulomb contribution 

Because we take different values for the nn and np 1 Sq scattering lengths, our calculations 
do include some higher-order isospin-breaking effects. Given that we have not done a consis- 
tent calculation to N 2 LO (or even NLO), the best we can say is that the ~160 keV difference 
between our result and that of Ref. 19j is comparable in size to the contributions of 70- 
110 keV that have been estimated using other approaches [3IH33J]. One caveat that should 
be made about these comparisons is that, as noted above, there is no model-independent 
way to separate Coulomb and other isospin-breaking effects in this EFT. 
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By adjusting the strength of the three-body force, we can fit the experimental binding 
energy. As in the triton case, we have done this for a range of cutoffs. The results are plotted 
in Fig. El They display the same limit-cycle behaviour as in the triton channel. Moreover 
the differences between the forces needed in the two channels are very small, reflecting the 
fact that a symmetric force can give a very good account of 3 He. 

To quantify the size of the isospin breaking in the three-body forces, it is better not to 
consider the values of H(A), since these cycle between ±00. The scale A or, equivalently A*, 
is a more appropriate measure. Its value of A = 382.46 MeV differs from that for the triton 
by 0.46%. The contribution of this force to the triton- 3 He splitting is about 50 keV. This is 



much larger than the 5 keV estimate of the three-body contribution in Ref. 33j] , which may 



reflect the fact that other isospin-breaking effects are missing from our calculatation. Also, 



the results of Ref. 



33( were obtained using an EFT that includes pion fields. 



V. SUMMARY 

We have studied 3 He within the framework of the pionless EFT, treating the Coulomb 
interaction nonperturbatively. We solve the set of integral equations developed by Kok et 



al. 20| , which can be thought of as extending the equation of Skornyakov and ter-M artirosian 
2l| to include the full off-shell Coulomb T-matrix. 

We find that a three-body interaction is required at leading order in 3 He, as in the triton. 
This force also exhibits the same limit-cycle behavior as a function of the cutoff A as found 



by Bedaque et al. js] for systems with purely short-range interactions. These results show 
that the Coulomb interaction is not singular enough to alter the deeply bound "Efimov" 
states in the three-nucleon system. 

We find that an isospin symmetric three-body force, fit to the triton binding energy, can 
give the 3 He binding to better than 1%. The scale parameter of the force that fits 3 He differs 
from the corresponding value for the triton by about 0.5%. This is of the expected size for 
an order-a electromagnetic effect. 

Our results demonstrate that the Coulomb interaction has no nonperturbative effect on 
the renormalisation of the three-body force in 3 He. Moreover the overall contributions of the 
Coulomb interaction to the binding energy of 3 He are small, perhaps unsurprisingly given 
the typical momenta involved, and so could have been calculated perturbatively. This will 
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not be the case if our approach is extended to describe low-energy pd scattering. However, 
as noted above, such an extension will require improved numerical techniques to handle the 
singularity of the Coulomb interaction. 

Our approach could also be used to determine the 3 He wave function from the correspond- 
ing eigenvector of the set of integral equations. This could then be used in calculations of 



electromagnetic properties of 3 He, based on the methods in Ref. 341 ] . 

Our present results do include certain other isospin-breaking effects beyond LO in this 
EFT, through the nn and np 1 Sq scattering lengths. A proper treatment of these in the three- 
nucleon system will require a consistent higher-order calculation. We have provided here 
the basic forms of the integral equations allowing for isospin-breaking in the NN scattering 
lengths. The relevant power counting for the terms needed in a full calculation can be 



determined using the techniques in Refs. 
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